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ABSTRACT 

The rotation curves of low surface brightness galaxies provide a unique data set 
with which to test alternative theories of gravitation over a large dynamic range in 
size, mass, surface density, and acceleration. Many clearly fail, including any in which 
the mass discrepancy appears at a particular length-scale. One hypothesis, MOND 
[Milgrom 1983, ApJ, 270, 371], is consistent with the data. Indeed, it accurately 
predicts the observed behavior. We find no evidence on any scale which clearly 
contradicts MOND, and a good deal which supports it. 
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1. 



Introduction 



There exists clear evidence for mass discrepancies in extragalactic systems. Application of 
the usual Newtonian dynamical equations to the observed luminous mass does not predict the 
observed motions. This leads to the inference of dynamically dominant amounts of dark matter. 

In a preceding paper we described the many difficulties that arise in trying to understand 
the data in terms of dark matter (McGaugh & de Blok 1998; hereafter paper I). In this paper 
we consider the alternative of a change in the fundamental equations of motion as the cause 
of the observed mass discrepancies. Many alternative theories of gravity have been posited to 
this end, but most have been ruled out (e.g., Sanders 1986). A few which may still appear to 
be viable are not in the light of data for Low Surface Brightness (LSB) galaxies. Length-scale 
dependent alteration of the inverse square law (e.g., the generic cases discussed by Liboff 1992 or 
the nonsymmetric gravity of Moffat & Sokolov 1996) can not explain the large variation in the 
scale on which the mass discrepancy appears (Fig. 3 of paper I) unless the length-scale of the 
theory is allowed to vary from galaxy to galaxy. Similarly, the linear potential theory explored 
by Mannheimm & Kazanas (1989) predicts rotation curves which should ultimately rise rather 
than remain flat. This is not obviously consistent with extant rotation curve data (Carlson & 
Lowenstein 1996), and the rotation curves of at least some LSB galaxies (e.g., UGC 128) remain 
flat well beyond the point where the upturn should be observed. Only the Modified Newtonian 
Dynamics (MOND) proposed by Milgrom (1983a,b,c) appears empirically viable, so we focus on 
testing it. 

The rotation curves of disk galaxies derived from gaseous tracers (Ha and H I) provide the 
strongest tests of alternative force-laws (e.g., Kent 1987; Begeman et al. 1991). With only the 
assumption of circular motion, it is possible to directly equate the centripetal acceleration 

V 2 

ac = li (1) 

with the gravitational acceleration 
determined from the Poisson equation 



L V «o / J 

(Bekenstein & Milgrom 1984). In no other type of system are tests so direct and free of 
assumptions. 
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(3) 



or any hypothesized alternative, for example 




(4) 



To test MOND we employ the data compiled in paper I (mostly from Broeils 1992, van der 
Hulst et al. 1993, and de Blok et al. 1996). These are rotation curves of disk galaxies spanning a 
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large range in size, luminosity, and surface brightness. Luminosity and surface brightness must 
trace mass and mass surface density modulo only the mass-to-light ratio of the stars in this 
alternative to dark matter. 

Previous tests of MOND have given mixed results. MOND does generally do a good job of 
fitting rotation curves (Begeman et al. 1991), but it is less clear that it works in systems other than 
disks (e.g., clusters of galaxies, The & White 1988). As emphasized by Milgrom (1983b, 1988), 
LSB galaxies provide particularly strong tests. Some attempts to test MOND with LSB dwarf 
galaxies found that it failed (e.g., Lake 1989). Milgrom (1991) pointed out some uncertainties and 
limitations in Lake's analysis, so it remains unclear whether these cases constitute contradictions 
of MOND. Other attempts to fit a limited number of LSB galaxies reported success (Milgrom & 
Braun 1988; Begeman et al. 1991; Sanders 1996). 

We test MOND with the substantial amount of new data we have accumulated on LSB 
galaxies (van der Hulst et al. 1993; de Blok et al. 1996). Section 2 describes the predictions 
of MOND relevant to LSB galaxies and tests each. Section 3 examines the viability of MOND 
in other systems. A summary of the results and some discussion of their implications is given 
in §4. Detailed MOND fits to the rotation curves of LSB galaxies are given by de Blok & 
McGaugh 1998 (paper III). Symbols and definitions follow the conventions of paper I. We adopt 
Hq = 75 kms -1 Mpc -1 throughout. 



2. The Modified Dynamics 

The Modified Newtonian Dynamics are an empirically motivated force law which can be 
interpreted either as a modification of gravity or of the law of inertia (see review by Milgrom 
1994). We are not interested here in which interpretation is preferable, but rather in testing 
whether MOND is indeed the correct force law. In this context it should be realized that MOND 
has not yet been developed into a full theory in the sense of General Relativity. However, as a 
force law it makes very precise and testable predictions just as the inverse square law did well 
before the adornment of subsequent theoretical elaboration. 

The empirical motivation for MOND is the observation that the rotation curves of high 
surface brightness (HSB) spiral galaxies are asymptotically flat. One must take care not to test 
only the fact the force law was constructed to realize. The properties of LSB galaxies were largely 
unknown at the time, and were not part of the input motivating MOND. Milgrom (1983b) stated 
that "disk galaxies with low surface brightness provide particularly strong tests" and made a series 
of specific predictions about LSB galaxies which constitute genuine tests of the MOND force law. 
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2.1. The Modified Force Law 

The MOND force law is 

a = \/gNao, (5) 

where is the usual Newtonian acceleration and oo is a universal constant. This applies only for 
a <C ao; for a > ao the behavior is purely Newtonian. There is an interpolation function /jl(x) with 
x = a/ao for connecting the two regimes (Milgrom 1983a). The function is required to have the 
asymptotic behavior n{x 3> 1) — ► 1 and fi(x <C 1) — > x. We will not specify //(x) here, and restrict 
our tests to the deep MOND regime where equation 5 is the effective force law. In practice this 
means x < 1 /2, where the precise form of fi(x) ceases to matter (i.e., [n(x) — x]/x < 10%). The 
MOND limit only occurs for extremely small accelerations: ao = 1.2 x 1CP 10 ms~ 2 (Milgrom Sz 
Braun 1988; Begeman et al. 1991). This is roughly 1 Angstrom per second per second, or 1CT 11 g®. 

For circular orbits about a point mass, Ai, in the MOND limit, 

v c 2 [gm 

This gives an asymptotically constant rotation velocity V c independent of R: 

V c A = a GM. (7) 

It is this behavior that gives rise to asymptotically flat rotation curves and the Tully-Fisher 
relation (Tully & Fisher 1977). 

Since MOND is a force law based on a modification at a particular acceleration-scale, the 
strongest tests of MOND are provided by systems with the lowest accelerations. The stars in HSB 
galaxies experience centripetal accelerations of order ao- Low surface brightness galaxies have very 
low mean accelerations, typically (a) w 10~ n ms~ 2 ~ ao/10. As a consequence, LSB galaxies 
provide strong tests of MOND. 



2.2. Testing MOND 

Solar system tests can be very precise, but do not probe the MOND regime. The transition 
between Newtonian and MOND regimes allows us to define a transition radius where a = ao, 

R t = GM/a . (8) 

For the sun, Rt ~ 0.1 light year. One would not expect to notice MOND effects until this radius, 
at which point the field of the Milky Way becomes significant. Classic solar system tests near 
to the sun are restricted to the regime a S> ao and only test the deviation of the interpolation 
function from its asymptotic limit (Milgrom 1983a). The same holds true for binary pulsars, which 
experience large accelerations well removed from the MOND regime. 



- 5 - 



Extragalactic systems are the only laboratories where the MOND regime is clearly probed. 
In fitting MOND to galaxy data, there are several adjustable parameters. One is the acceleration 
constant itself. This must be a universal constant; once measured accurately in one galaxy the same 
value must apply in all galaxies. We adopt ao = 1.2 As -2 and keep it fixed. This was determined 
by Begeman et al. (1991) to be the best fit to a sample of galaxies with high quality H I rotation 
curves which has no overlap with our new LSB galaxy data. Slight adjustments to the distance 
to a galaxy can sometimes improve a MOND fit (Begeman et al. 1991). This occurs because of 
the interplay between the MOND mass, which depends only on the distance-independent circular 
velocity, and the gas mass, which varies as D 2 . We keep D fixed at the value determined from the 
assumed Hq. For other choices of Hq there is some limited freedom to adjust ao to compensate. 
As with ao and D, adjusting the inclination of a disk galaxy can improve some fits. Inclination is 
crucial since MOND masses depend on V 4 / sin 4 (i). In this paper we keep i fixed; see paper III for 
further examination of this point. 

So far we have named three parameters which could in principle be adjusted. In practice, 
there is very little freedom to do this, and we will keep all of them fixed. The only truly free 
parameter is the conversion factor from light to mass, Y*. We will treat this as a value to be 
determined and compared to the expectations for stellar populations. 

In the following subsections, we discuss particular tests of MOND with LSB galaxies. 



2.3. The Tully-Fisher Relation 

A very strong prediction of MOND is a single universal Tully-Fisher relation. 11 The relation 
between asymptotic velocity and the mass of the galaxy is an absolute one." (Milgrom 1983b). 
This follows from equation (7) which gives 

V c 4 = a GM = a G (^^j L oc T b L (9) 

where the baryonic mass-to-light ratio Y& = Mb/L. Note that R does not appear in this equation, 
nor do the dark matter galaxy formation parameters of equation (13) in paper I. 

That R does not appear in equation (9) is of fundamental significance. In the case of purely 
Newtonian dynamics, V 2 = GAi/R. Since LSB galaxies have, by definition, larger R than HSB 
galaxies of the same luminosity, they should not fall on the same Tully-Fisher relation. Yet they 
do (Sprayberry et al. 1995; Zwaan et al. 1995). Invoking a length-scale for the mass which is 
different from that of the light helps not at all: serious fine-tuning problems occur in the dark 
matter picture (paper I). In MOND the Tully-Fisher relation follows simply and directly from the 
form of the force law. 
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The error budget for scatter in the Tully-Fisher relation in magnitudes is 



\5M\ 2 = 18.86 



sv c 


2 




+ 1.18 


v c 





Comparison of this with equation 12 of paper I shows that the term involving the central surface 
brightness has disappeared. The real difference is larger, though. In the dark matter picture, 
there should be some intrinsic variation in V c from halo to halo. In MOND, mass and velocity 
are strictly related; the term for the error in V c in equation (10) now refers only to observational 
uncertainties and not also to intrinsic variations. All but one of the terms which should contribute 
intrinsic variance in the dark matter picture disappear. To obtain a single universal Tully-Fisher 
relation with little scatter all that is really required in MOND is that Mb/L be roughly constant: 
SM = 1.086 |(5Tft/Tb|. Small intrinsic scatter in the Tully-Fisher relation is much easier to 
understand with MOND than with dark matter. 

The slope of the Tully-Fisher relation, L oc V y with y = 4 (Aaronson et al. 1979; Tully & 
Verheijen 1997; paper I), is dictated by the form of the force law in MOND. No particular slope 
is required in the case of dark matter. A slope of y = 4 is sometimes attributed to the virial 
theorem (e.g., Silk 1997), but a slope closer to y = 3 has also been suggested as the value arising 
in plausible disk galaxy formation scenarios (e.g., Mo et al. 1997). In any case, a constant disk 
plus halo mass surface density must be arranged. Multi-component disk-halo models generally 
predict that some signature of the disk component should be visible in Tully-Fisher residuals, but 
none are present (Courteau &; Rix 1997). More generally, there is no signature of the transition 
from disk to halo domination (the "disk-halo conspiracy" ) . These all occur naturally if there is no 
halo and the disk is the only component, as in MOND. 

The Tully-Fisher relation of LSB galaxies provides a new test of MOND. That HSB galaxies 
obey the Tully-Fisher relation was known when MOND was developed, so this does not itself 
constitute a test. That LSB galaxies fall on the same relation with the same normalization was 
not known. That they should constitutes a genuine prediction (Milgrom 1983b): "We predict, for 
example, that the proportionality factor in the M. oc V£ relation for these galaxies is the same as 
for high surface density galaxies." 



2.4. The T-S Conspiracy 

If the MOND force law is the basis of the Tully-Fisher relation, it should be possible to derive 
from MOND the T-S conspiracy inferred with conventional dynamics (Zwaan et al. 1995; paper 
I). This requires T Q 2 S oc constant, where T Q is the conventional dynamical mass-to-light ratio. In 
MOND, the typical acceleration of a disk is proportional to the square root of its characteristic 
surface density so that 

(a) 2 oc a = T 6 S (11) 
(equation 6 of Milgrom 1983b and equation 19). 
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The severity of the inferred Newtonian mass discrepancy M^/Mb depends on the extent to 
which the acceleration is in the MOND regime: 

M N 1 Y G a 



The approximation effectively becomes an equality for (a) < clq/2. As a result, T G = Tj,ao/(a). 
Using equation (11), this becomes Y c 2 oc (ag/c^Y 2 , = agTb/S, i.e., 



Another test of MOND is provided by the values required for the mass-to-light ratio of stellar 
populations, Y*. These should be consistent with what we know about stars. Since there is no 
dark matter, the kinematic data provide a direct measure of the luminous mass. 

The only free parameter in equation (9) is Y*. The constants ao and G are known, V c , L, 
and the gas mass M. g are measured. Hence, = M. — M. g and Y* = Ai*/L where M. is the 
total dynamical mass indicated by MOND. For spiral disks we expect Y* in the range of 1 - 2 in 
the .B-band, with the range 0.5 < Y* < 4 being credible (e.g., Larson & Tinsley 1978; Bruzual & 
Chariot 1993). 

The gas mass is inferred directly from the 21 cm flux multiplied by a conversion factor to 
account for helium and metals: M g = IAMhi- The conversion factor is probably closer to the 
primordial value 1.32 for LSB galaxies since these systems are very metal poor (McGaugh 1994; 
Ronnback & Bergvall 1995). This difference is very small compared to the accuracy with which 
Y* can be estimated. We assume molecular gas is not a major mass component. High surface 
brightness galaxies are gas poor with molecular gas mass being a small fraction of the stellar 
mass (Young & Knezek 1989; McGaugh & de Blok 1997). Low surface brightness galaxies are H I 
rich but apparently have very little molecular gas (Schombert et al. 1990; de Blok 1997). Any 
molecular mass which is present will be attributed to stars, resulting in Y* which are slightly too 
high. Again, this is a small effect compared to the uncertainty in Y*. 

Strictly speaking, equation (7) holds precisely only for point masses. Disk galaxies are 
reasonably well approximated by an exponential mass distribution for which 



M b v{x) T b (a)' 



(12) 



Y Q 2 S oc Y^Oq ~ constant. 



(13) 



The Y-E conspiracy does follow from MOND. 



2.5. Stellar Mass-to-Light Ratios 




(14) 



and 



g N (R) = ttGcjo- [I (R/2h)K Q (R/2h) - / 1 (i?/2/ l ) J PC 1 (i?/2/t)] , 



(15) 
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where /„ and K n are modified Bessel functions of the first and second kind (Freeman 1970). 
Combining these and following the MOND prescription a = ^/qnOo gives a dimensionless factor 
x(R) correcting the simple point mass formula for the disk geometry: 

2h 3 \l-e- R / h (l + R/h)} 
X(R) = ^ rrfnl nu^ fnl nu. r——±———i (16) 



R 3 [I {R/2h)K (R/2h) - I 1 (R/2h)K 1 (R/2h)} 
so that 

M(R)= X (R)^ 1 - (17) 

This is a fairly mild correction, as x(R) does not deviate far from unity at most radii. For 
consistency with paper I we evaluate it at R = Ah, where x(4/i) = 0.76. The effect is to reduce the 
mass somewhat as a disk geometry rotates faster than the equivalent spherical mass distribution. 

The stellar mass is computed as 

V A 

M* = 0.76^- - IAMhi- (18) 

We have determined Y* in this fashion for all galaxies meeting the requirements stipulated in 
paper I (essentially everything for which we have a measurement of V c , L, jjlq, and h). The results 
are given in Table 1 and plotted in Fig. 1, which is analogous to Fig. 4 of paper I. 

The stellar mass-to-light ratios inferred from MOND are consistent with those expected for 
disk population stars. In the mean, (Y*) = 1.9, and T* = 1.6 in the median. There is a fair 
amount of scatter about the mean value. No trend of Y* with no is obvious in Fig. 1(b). Slight 
trends of Y* with Mb and h are perceptible in Figs. 1(a) and (c), but are not highly significant as 
they are strongly influenced by the low Y* values of a few small galaxies. The directions of these 
trends are nevertheless consistent with those of color-magnitude and color-size relations. 

MOND can be quite sensitive to the parameters we have held fixed: do, D, and i. It has 
not been necessary to adjust any of these parameters to get reasonable results from the global 
approach taken here. By global, we mean that we use only the asymptotic flat circular velocity V c 
to determine M. In paper III we present fits to the full shape of V(R). This is a different exercise 
which gives slightly different results for the mass-to-light ratios of individual galaxies. 

The precise distance to a galaxy is sometimes important because of the different -D-dependence 
of gas mass and total MOND mass. Total mass M. oc is measured independent of distance 
(for small redshifts), but M g oc D 2 . An overestimate of the distance can lead to a situation in 
which it appears that the gas mass exceeds the total dynamical mass. For a sample with many 
gas rich galaxies, it would be surprising if this did not occasionally manifest itself. Indeed, it 
is well known in the case of DDO 154 (Milgrom & Braun 1988; Begeman et al. 1991; Sanders 
1996). Examination of Table 1 shows that F571-V1 and F574-1 might be similar cases, as both 
are inferred to have unreasonably small or negative Y*. A 20% reduction in the assumed distance 



-9- 



to these galaxies would have the effect T*(F571-V1) -> 0.5 and Y*(F574-1) -> 0.8. Given the 
uncertainty in the distance scale, neither case provides a clear counterexample to MOND. 

MOND is also very sensitive to the inclination since the observed line of sight velocity must 
be corrected by sin~ 4 (i). Inclination determinations can be difficult and are sometimes quite 
uncertain for LSB galaxies (McGaugh & Bothun 1994; de Blok et al. 1995; de Blok et al. 1996). 
In the case of F571-V1, Y* becomes positive if i = 35 — > 34. We certainly would not claim to be 
able to determine inclinations this accurately by information independent of MOND (paper III). 

In the case of F574-1, it is not obvious that the asymptotically flat part of the rotation curve 
has been reached (de Blok et al. 1996). If V c is underestimated, the MOND mass will be too small. 
The stellar mass to light ratio becomes more reasonable (T* > 0.5) if V c = 91 — > 99 kms -1 , which 
is possible within the errors. It is also possible that the factor x(i?) we have assumed could be an 
underestimate in this case. 

In sum, MOND yields reasonable values for the mass-to-light ratios of stellar populations. 
This is not a trivial result. Since Y* oc V^, even small errors in V c can make Y* seem incorrect. 
There are a few individual galaxies for which Y* may seem a bit unreasonable, but is this 
surprising given the nature of astronomical data? 



Since there is no dark matter in the MOND hypothesis, mass surface density must be 
correlated with luminous surface density. A way to see this without invoking stellar populations is 
through the characteristic parameter £ defined by Milgrom (1983b) for exponential disks: 



This is proportional to the characteristic acceleration V 2 /h of a disk and to the square root of its 
mass surface density (a oc M/h 2 ). 

By using the observed velocity and scale length, we can construct the quantity £ which 
provides a direct dynamical estimate of the mass surface density. Values computed for £ are given 
in Table 1. The mass surface density indicated dynamically by £ must be correlated with with 
the mass surface density of ordinary matter, indicated by the central surface brightness jiQ. These 
quantities are plotted against each other in Fig. 2. There is a strong correlation, with regression 
coefficient 1Z = 0.82. The slope and normalization of the relation is consistent with that expected 
in MOND. From equation (19), one can derive 



assuming an exponential mass distribution with M = 2-iraoh 2 . A formal fit to the data gives a 
slope of —0.21, consistent with the expected —0.2. The normalization is also correct for Y& « 4, a 



2.6. 



Mass Surface Densities 




(19) 




(20) 
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reasonable value given the assumption of an exponential mass distribution. Deviations from this 
idealized case will alter the factor 2ir in equation (20) and hence the precise value of Y& that is 
inferred. 

2.7. Conventional Quantities 

The effects of MOND, when interpreted in terms of conventional dynamics, lead to specific 
predictions for the behavior of conventional quantities. For example, the Newtonian mass-todight 
ratio in any given galaxy should show no mass discrepancy until beyond the transition radius 
defined by equation (8). After this point, T a (R) should increase as R increases. As one examines 
galaxies of decreasing central surface brightness, two effects should be apparent: the severity of 
the discrepancy should grow, and the radius (measured in scale lengths) at which the discrepancy 
sets in should shrink. If at every point in a galaxy a(R) <C do, the mass discrepancy should be 
apparent at essentially R = (Milgrom 1983b). 

In Fig. 3 we plot T G against R/h for all available galaxies (Tables 1 and 2 of paper I). 
The mass-todight ratio is computed as the accumulated mass within radius R divided by the 
accumulated luminosity at the same point. Mass is assumed to be distributed spherically 
[M oc RV 2 (R)] which is an adequate approximation for present purposes if not strictly true at 
small radii. We have subdivided the data in Fig. 3 into one-magnitude bins in surface brightness. 
The rise of T D does appear to become more pronounced as no declines. There is no clearly defined 
transition radius, but all disk galaxies have a <; do- This may be a requirement for disk stability 
(see §3.3). 

What really matters is the acceleration. Another way of stating the prediction of MOND is 
that (a) should decrease with surface density (equation 19). In Fig. 4 we plot the rotation curves 
as the centripetal acceleration required to produce them: a = V 2 (R)/R. There is some scatter 
induced by inclination uncertainties and intrinsic scatter in Y*. Nevertheless, there is a clear trend 
for (a) to decline with surface brightness, as expected in MOND. 

2.8. Rotation Curve Shapes 

The shapes of the rotation curves of exponential disks is controlled in MOND by the 
parameter £ which is closely related to the surface brightness. Milgrom (1983b) predicts "a 
correlation between the value of the average surface density (or brightness) of a galaxy and the 
steepness with which the rotational velocity rises to its asymptotic value. Small surface densities 
imply slow rise of V." We noted in paper I the slow rate of rise of the rotation curves of LSB 
galaxies (fact 2). Our empirical statement is essentially identical to the prediction of Milgrom 
(1983b). As a measure of this effect, Milgrom (1983b) suggested the radius (in scale lengths) at 
which V(R/h) = V c /2. This contains the same information as R34 in paper I which should thus be 
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correlated with surface brightness. That is, Fig. 12(b) of paper I is in fact a prediction of MOND. 

We illustrate this further by plotting the rotation curves logarithmically to show their shape 
(Fig. 5), binned by surface brightness as before. The rotation curves of the highest surface 
brightness galaxies rise rapidly, often reaching the asymptotic velocity within 1 scale length. The 
rate of rise becomes slower as surface brightness decreases. The asymptotic flat velocity is only 
reached at 2 or more scale lengths in the lowest surface brightness galaxies, as expected in MOND. 

Another systematic of rotation curves is that the rate of rise of V(R/h) is well correlated with 
absolute magnitude (Fig 12a of paper I). That this occurs is not predicted or required by MOND. 
Hence it neither supports nor contradicts MOND. 

Though the rotation curve shapes are consistent with MOND, the predictions discussed in 
this section are less strong than previous ones. The systematic trend of -R34//1 with £ is predicted 
to be fairly weak, and the presence of a bulge component can have an additional effect on the 
shape of V(R) not accounted for by the quantity £ defined for exponential disks. Though bulges 
are generally anemic in LSB galaxies, they are present in the higher surface brightness galaxies in 
Fig. 5. The question then becomes whether MOND can fit the rotation curves in detail given the 
actual observed luminous mass distribution and not just the exponential disk approximation. 

2.9. Residuals of Rotation Curve Fits 

MOND is known to work well in fitting rotation curves, particularly in HSB galaxies 
(Begeman et al. 1991; Sanders 1996). Except for a few dwarfs in those samples, V(R) reaches the 
asymptotic flat velocity very rapidly. This is not true in LSB galaxies: MOND must not only give 
an asymptotically flat rotation curve, but it must also give the gradual observed rise. It was not 
contrived to do this, so detailed fits to LSB galaxy rotation curves provide strong tests. 

An example of a detailed MOND fit to an LSB galaxy is given in Fig. 6. Fits for the entire 
sample are the subject of paper III. We make use of those results here with only a few comments. 
The prediction of V{R) follows with one fittable parameter, Y*. The freedom in adjusting Y* 
is very limited. It does not have a strong effect on the fit (V oc Y* 1 / 4 ) and it must return a 
reasonable value for a stellar population. The shape of the stellar mass distribution is constrained 
to be the same as that of the light; T* acts only as a normalization factor and is not allowed 
to vary with radius. The gas mass is comparable to the stellar mass in these galaxies. There 
is very little uncertainty in the conversion factor from 21 cm flux to H I mass, and we have no 
freedom to adjust it. The importance of the gas component often constrains Y* so tightly that in 
many cases the "fits" are effectively parameter-free. For nine LSB galaxies, good fits were found 
immediately with one parameter (Y* only) fits. For six others, achieving tolerable fits required 
adjustment of the inclination (paper III). MOND is very sensitive to this since it enters through 
sin~ 4 (i) and a number of the LSB galaxies in our sample are fairly face-on. Two parameter (T*, 
i) MOND fits remain tightly constrained, and both parameters are subject to independent checks. 
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In contrast, dark matter fits require at least three parameters (T* and two halo parameters). The 
two halo parameters are not constrained by any data independent of the rotation curves, leading 
to notorious degeneracies (Athanassoula et al. 1987; Kent 1987; de Blok & McGaugh 1997). 

The residuals to the MOND fits to the rotation curves of LSB galaxies from paper III are 
shown in Fig. 7 as a function of the critical parameter x = a/aQ. The data extend over an order 
of magnitude range, all with x < 1 and most with x < 1 /2- (See Sanders 1996 for fits to galaxies 
with x > 1 as well as x < 1). Each point is a measured point from a rotation curve; the residuals 
of all LSB galaxy fits are plotted together. The data closely follow the line of zero residual with 
~ 10% scatter, about what is expected for our errors. In absolute terms, the mean deviation 
AV = V obs - Vmond = 0.6 ± 5.6 kms" 1 . 

The strongest test is provided by the points with x < V2 where the form of the interpolation 
function fj,(x) is insignificant. In this part of the diagram, not only is the mean residual small, but 
there is no systematic trend of the residuals with x. For x > 1 /2, there is a hint of a systematic 
deviation. The difference is not highly significant, but it might indicate that the actual from of the 
interpolation function is slightly different than the one assumed. The explanation may be more 
prosaic, as some of the data are modestly affected by beam smearing (see de Blok & McGaugh 
1997 for an extensive discussion). 

There is one further consequence of Fig. 7. Even though the mass discrepancy does not 
appear at a particular length scale, it does appear at a particular acceleration scale. MOND is the 
effective force law in disk galaxies. 

3. Other Evidence 

The data for LSB galaxies are consistent with MOND. Indeed, each of the specific predictions 
Milgrom (1983b) made about LSB galaxies is confirmed. This means something. Is it possible 
that MOND, not dark matter, is the solution to the mass discrepancy problem? 

There are many other systems besides disks in which the mass discrepancy is apparent. 
If MOND is correct, it must also be consistent with tests from these other systems. Tests are 
provided by any system with (a) <C ao; in systems where (a) S> ao, no mass discrepancy should be 
inferred unless there is genuine hidden mass. Here we review the evidence provided by other types 
of systems on different scales. 

3.1. Dwarf Spheroidals 

A class of galaxies which are low in surface brightness but quite distinct from the LSB galaxies 
we have so far discussed are the dwarf Spheroidals (dSph) found in the Local Group. They 
are distinct for a number of reasons, most importantly that they are three dimensional systems 



-13- 



supported by the velocity dispersion of their constituent stars, not rotating disks. As LSB systems, 
they should have low accelerations and provide a strong test of MOND: "Effects of the modified 
dynamics are predicted to be particularly strong in dwarf elliptical galaxies" (Milgrom 1983b). 

That dSph galaxies have large mass discrepancies is clear (Mateo et al. 1993; Vogt et al. 1995; 
Mateo 1996). However, it is less clear that the discrepancies are of the correct magnitude for 
consistency with MOND. Gerhard & Spergel (1992) and Gerhard (1994) have applied this test. In 
some cases, MOND seems to do well, giving stellar mass-to-light ratios which are reasonable for 
the stellar populations of these galaxies. In other cases, the results are less satisfactory. 

Dwarf Spheroidals provide a strong test in the sense that the acceleration is low, comparable 
to LSB galaxies. However, they do not provide as clean a test. What one really needs to know is 
the three dimensional velocity ellipsoid. All that is possible, of course, is to obtain a line-of-sight 
velocity dispersion <j. Additional assumptions are therefore required to interpret the observations. 
We must assume that the MOND equivalent of the virial theorem applies (Milgrom 1984; Gerhard 
& Spergel 1992; Milgrom 1995) and that each system is either spherical or isotropic. 

In general, the mass of an isolated system in the MOND regime is 

9(V 2 ) 2 

M = - A ^%- (21) 
4 a G 

(Milgrom 1995). The assumption of sphericity or isotropy is necessary to relate the rms velocity 
to the observable line-of-sight velocity (V 2 ms ) = 3? 2 . Thus the effective equation is 

M = (22) 
4 a G 

Note that the assumed relation between rms and line-of-sight velocity has a significant impact 
in the geometric constant relating mass to velocity. Since A4 oc ^ 4 , a 20% error in <; will lead to 
a factor of two error in mass. This need not arise as an error. It might occur even with perfect 
observations simply because the true velocity ellipsoid may not have the form we are obliged to 
assume. 

Gerhard & Spergel (1992) and Gerhard (1994) concluded that MOND fails because it yields 
unreasonable Y* for some dwarf Spheroidals, especially Fornax and Ursa Minor. By reasonable, 
we mean 1 < Y* < 6 in the F-band. There is considerable uncertainty in the stellar content of 
these galaxies. Originally thought to be old implying high Y*, at least some appear surprisingly 
young (e.g., Mateo et al. 1991; Smecker-Hane et al. 1994). The T* of Gerhard (1994) are plotted 
in Fig. 8. One can immediately see two things: there are no error estimates, and many of the 
values are in fact quite reasonable. So, do the unreasonable values falsify MOND, or are the errors 
simply very large in those cases? 

Data of the sort required for this analysis have accumulated rapidly in the period during 
which these analyses were performed. As a result, rather better data now exist than did at the 
time of these analyses. Milgrom (1995) has revisited this issue. The results he obtained using the 
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same method and equations as Gerhard & Spergel (1992) are also plotted in Fig. 8. Within the 
errors, the results give reasonable values for T*. The only thing that has really changed is that 
the data have improved — data which come from the same sources. 

The data continue to improve, so we repeat this analysis using the compilation of data 
provided by Mateo (1996). Equation (22) above only applies for isolated spheres. When a system 
is subject to a dominant external field (e.g., that of the Milky Way), the MOND behavior is 
quasi-Newtonian 

M = 2G-j f c, 2 R v (23) 

with a modified effective constant of gravitation, G e ff — > aoG/(a). Here, Ry is the virial radius 
(equation 10 of Gerhard & Spergel 1992). Which equation applies depends on the relative 
accelerations imposed by the internal and external fields: 



3<? 2 /2ifc 
Kw/D 



V = 9^ (24) 



where Rc is the core radius of the King profile which fits the light distribution of the dSph in 
question, D is its galactocentric distance, and Vmw is the asymptotic velocity of the Milky Way. 
For the latter we adopt Vmw = 220 kms -1 . If rf > 1 the internal field dominates and the object 
can be treated as isolated (equation 22). If r\ < 1, the external field due to the Milky Way is 
dominant, and the behavior is quasi-Newtonian (equation 23). As can be seen in Table 2, most 
dSph galaxies for which good kinematic data exist are probably influenced by the external field of 
the Milky Way. When 77 « 1, neither of these limits are really appropriate. This is the case for 
many objects when the uncertainties in the many parameters determining 77 are considered. A 
significant external field can elongate the internal field and hence affect the assumed geometry, so 
there are significant uncertainties besetting this analysis regardless of the precision of the data. 
Only Leo I provides a test in a reasonably isolated system. 

Where Mateo (1996) gives different values for we adopt the one based on the greater 
number of stars. In most cases, <r has only been measured in the central regions. A global mean is 
more appropriate, so in the case of Fornax where a velocity dispersion profile is available we use it 
to estimate s « 12 kms -1 . We do not compute T* for two dwarfs for which velocity dispersions 
are listed, LGS 3 and Sagittarius. The velocity dispersion listed for LGS 3 gives Y* = 2.5 but is 
based on only 4 stars, so the uncertainties render this meaningless. In the case of Sagittarius, 
the very close interaction with the Milky Way undermines all of the assumptions on which the 
analysis is based. 

Results of our analysis of the data compiled by Mateo (1996) are listed in Table 2 and shown 
in Fig. 8 together with the previous determinations of Gerhard (1994) and of Milgrom (1996). 
Our error bars are based only on the error in <j given by Mateo (1996); they do not include the 
large uncertainties in the luminosity, surface brightness, and distance to these galaxies. Neither do 
we attempt to estimate any of the systematic uncertainties discussed above. Hence the error bars 
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plotted in Fig. 8 are a lower limit to the true range of uncertainty. Nevertheless, they do show 
that our determinations are consistent with previous ones, when the errors are considered. 

Most dwarf Spheroidals do in fact lie in the range 1 < Y* < 6 according to all three 
independent analyses. Marginal cases are Draco and Ursa Minor (Y* too high) and Fornax (Y* 
too low). The deviation of Draco and Fornax from the reasonable range is marginal (~ la). 
Draco has rj very near to unity which will tend to cause a modest overestimate of Y* (Milgrom 
1995). Ursa Minor is the most problematic case, but even there the deviation is only ~ 2a. There 
is a complication in this case in that some of the velocity dispersion may result from rotation 
(Mateo 1996). This can affect the analysis since rotational velocity enters the calculation with 
less weight than does velocity dispersion. Also, Mateo (1996) lists two somewhat different velocity 
dispersions for this object, both based on a sizable number of stars. Armandroff et al. (1995) give 
c; = 8.8 ± 0.8 kms -1 while Hargreaves et al. (1994) give <; = 6.7 ± 1 kms -1 . Using the latter value 
gives Y*(Ursa Minor) = lOj^, not significantly unreasonable. 

Another complication is the systematic uncertainty in the parameters of the Milky Way. 
Since rj is close to unity in most cases, a difference in the assumed strength of the external field 
of the Milky Way might make a difference. It has been suggested that Vmw could be as low as 
180 kms^ 1 (Oiling 1997, private communication). If we re-do the analysis with this value, little 
changes. In most cases there is a slight shift towards lower Y*. The exception is Ursa Minor, 
where it makes a substantial difference. For the higher velocity dispersion of Armandroff et al. 
(1995), Y* = 17 — ► 10. In principle, uncertainties in the distance to individual dwarfs can have 
similar effects. 

Gerhard (1994) stresses that Y* varies by a factor of 20 between Fornax and Ursa Minor. It 
is true that this is a lot, but it is also true that these are the two most extreme points on either 
side of a sensible mean value. This is bound to happen at some level, and the degree to which it 
occurs here does not seem outrageous given the uncertainty in the data and the analysis. 

In sum, there is no evidence which clearly contradicts MOND in the data for dSph galaxies. 
Agreement between these data and MOND has generally improved as the data have improved. 
The qualitative prediction of MOND, that dSph galaxies should have high mass discrepancies, is 
not in dispute. MOND also appears to work quantitatively. 

3.2. Giant Elliptical Galaxies 

The kinematic data which exist for giant elliptical galaxies is rather limited. This is because 
observations of these systems are generally restricted to absorption lines which can typically 
be measured only within one half light radius (R e ) where the surface brightness is high (e.g., 
Bertin et al. 1994). Application of these data as a test of any hypothesized force law is limited 
by the limited range of radii probed. It is further hindered by the fact that ellipticals are quite 
complicated systems. The brightest ellipticals are predominantly pressure supported, but often 
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have significant rotation as well (Davies et al. 1983). Even in the absence of rotation, anisotropics 
in the velocity ellipsoid must exist to give rise to the observed range of shapes of ellipticals. This 
is quite problematic for applying equation (22): we simply do not have adequate knowledge of the 
geometry of the orbits. 

For testing MOND, an even more severe constraint applies. The typical acceleration at radii 
for which measurements are available (Bertin et al. 1994) is 

3 <? 2 

a ~2«~^ 3a0 ' ^ 
Extant data do not probe the MOND regime. 

Giant ellipticals therefore provide no direct test of MOND. They do provide several indirect 
tests, however. One is that there should be no apparent mass discrepancy where a > ao- This is 
consistent with the persistent lack of a clear need for dark matter in elliptical galaxies (van der 
Marel 1991; Bertin et al. 1994). 

Though the data may not probe the MOND regime, MOND effects must matter at some level 
in elliptical galaxies. One might expect some critical phenomenon associated with the MOND 
scale ao- In the simple case of spherical galaxies in the MOND limit, Milgrom (1983b) expects a 
Faber- Jackson (Faber & Jackson 1976) relation of the form Ai oc ? 4 . If ellipticals are approximately 
isothermal, then Milgrom (1984) also expects a Fish (1964) law. MOND may therefore manifest 
itself in the regularity of the Fundamental Plane (Djorgovski & Davis 1987). The analog of the 
Fundamental Plane for disks is the Tully-Fisher relation. With //q or Re as a third parameter, 
the "fundamental plane" of disks is perpendicular to the luminosity-velocity plane, resulting in a 
narrow Tully-Fisher relation for galaxies of all surface brightnesses. Similarly, the Fundamental 
Plane of ellipticals is viewed nearly edge-on in the luminosity-velocity dispersion plane. This 
suggests that the same effect is at work. That there is a modest tilt to the Fundamental Plane 
merely indicates some modest systematic trend of Y* or the degree of pressure support with 
luminosity. 

A generalization of the Fundamental Plane by Burstein et al. (1997) finds that essentially all 
extragalactic systems form narrow structures in "«-space." This is indicative of a universal scale 
like ao- The K-structures exhibit little surface density dependence, the signature of the MOND 
force law. 

Elliptical galaxies with shell systems provide an interesting probe to large radii. Hernquist & 
Quinn (1987) use the remarkably regular shell system NGC 3923 to argue that MOND predicts 
the wrong number of shells for the periods implied by the observed relative radii, and also that 
the slope of the radius-shell number relation is incorrect (their Fig. 2). The latter argument is 
based on normalization to the outermost shell (Ngf^n = 1). This is a little unfair, as only the 
shape is predicted. The predicted shape is approximately correct for 1 < N s h e ii < 17 — it is 
the outer shell which is deviant. Unlike the majority of the shells in this system, this outermost 
shell is substantially offset from the major axis. Alignment with the major axis is central to the 
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argument, which assumes the shell system was created by a simple radial phase-wrapping merger 
event (Quinn 1984). Any deviations from this simple idealized case complicates the interpretation. 

The periodicity argument Hernquist & Quinn (1987) make is more persuasive. For a phase 
wrapped system, MOND should have formed more shells than they count. This requires that 
phase wrapping is an adequate model and that nearly all shells have been detected. Hernquist 
& Quinn (1987) argue that it would be difficult to miss the additional shells that seem to be 
required by their MOND argument. However, Prieur (1988) did discover additional shells which 
had previously been missed. 

Though a good argument, the period-number relation of Hernquist & Quinn (1987) does not 
constitute the direct test of the MOND force law. The argument is based on a very idealized 
realization of simple phase- wrapping as the result of a minor merger. Prieur (1988) shows that 
NGC 3923 is not as clean a system as required by this picture. Indeed, Hernquist & Spergel 
(1992) suggest that such rippled systems can result from major mergers (as already suggested 
from observational evidence by Schweizer 1982 and McGaugh & Bothun 1990). They specifically 
cite NGC 3923 as a likely example. Even in the conventional context, they abandon the minor 
merger phase wrapping hypothesis which is the basis of the argument against MOND. 

Another interesting argument involving an elliptical galaxy is the apparent difference between 
the shape and orientation of the optical and X-ray isophotes in NGC 720 (Buote & Canizares 
1994). These authors point out that in any modified theory of gravity, the isopotential surface 
presumably traced by the gas should not differ from that determined by the dominant stars. Any 
significant difference would impose a geometrical requirement for a dark mass component. 

In other elliptical galaxies Buote & Canizares (1994) examine there is no apparent difference 
in this sense. However, in the case of NGC 720, the stars have an ellipticity e* = 0.4 and the gas 
e g = 0.2 - 0.3 with a difference in position angle of 30 ± 15° (Buote & Canizares 1994, 1996). The 
X-ray isophotes twist and appear to become more pointy than the optical isophotes outside ~ lR e - 
Buote & Canizares (1994) argue that the potential due to the stars can only become rounder with 
increasing radius, so the pointy X-ray isophotes provide geometrical evidence in favor of dark 
matter. 

This situation is rather puzzling even with dark matter, since the stars should still contribute 
substantially to the mass over the observed region. Even though the X-ray isophotes appear to 
become more elongated than those of the stars, Buote & Canizares (1997) derive a dark matter 
mass potential which is more round than either stars or gas. This is less easy to accomplish in 
MOND since there is no freedom to vary the shape and position angle of the dominant mass as 
there is with dark matter. Nevertheless, both dark matter and MOND imply potentials consistent 
with or slightly rounder than the isophotes of the stars and rounder than the X-ray isophotes. The 
interpretation of NGC 720 seems difficult with either dark matter or MOND if we accept the data 
at face value. 

The basic argument of Buote & Canizares (1994) is based on the difference between the X-ray 
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and optical isophotes. The X-ray isophotes are very ragged (see Fig. 2 of Buote & Canizares 1994 
or Fig. 1 of Buote & Canizares 1996) and do not obviously provide a strong constraint. Paper III 
provides several examples of how misleading the shapes of ragged isophotes can be. 

The geometrical argument of Buote & Canizares (1994) is valid in principle. Observing many 
more galaxies with considerably higher signal to noise would prove interesting. At the present 
time, these and other data for giant elliptical galaxies provide only weak and indirect tests of 
MOND. None of these data clearly contradict it. 

3.3. Disk Stability 

Another indication of a mass discrepancy is in the long term existence of dynamically cold 
spiral disks. Purely Newtonian disks are subject to global instabilities which rapidly lead to their 
demise unless stabilized by a dominant dark halo (Ostriker & Peebles 1973). What really matters 
here is the ratio of binding to kinetic energy; this can be achieved either with dark matter of by 
altering the force-law. Milgrom (1989) showed analytically that MOND disks are somewhat more 
stable than purely Newtonian ones. 

It turns out to be very difficult to adopt standard N-body codes to address this problem 
properly (Mihos, private communication). Brada (1996) developed an alternative approach based 
on the multigrid algorithm. This supports and extends the analytic conclusions reached by 
Milgrom (1989): MOND disks are more stable than purely Newtonian disks both locally and 
globally. The additional stability is fairly modest, roughly equivalent to that provided by a halo 
with M(R = 5h) « 3M disk (Brada 1996). 

This leads to a further test. The stability properties predicted by dark matter and MOND 
diverge as surface brightness decreases. In the case of dark matter, the halo mass enclosed by the 
disk increases systematically with decreasing surface brightness (Fig. 4b of paper I). Self-gravity in 
the dynamically cold disk is the driving force for bars and spiral arms, but becomes progressively 
less important. At some point such features should be completely suppressed by the dominant hot 
halo. In contrast, the amount of additional stability provided by MOND depends only weakly on 
surface density and the self-gravity of the disk always matters. 

Spiral features appear feeble in LSB galaxies, but are clearly present (Schombert et al. 1992; 
McGaugh et al. 1995; de Blok et al. 1995; Impey et al. 1997). Brada (1996) predicts the growth 
rate of the m = 2 mode for both dark matter and MOND (his Fig. 11). Using the observed 
accelerations as a scale, LSB galaxies correspond to his mass models with M. < 0.1 Below this 
level, instabilities are almost completely suppressed in the dark matter case (see also Mihos et al. 
1997). 

In order to generate spiral structure internally, the disk needs to be rather heavy (Athanassoula 
et al. 1987). In LSB disks, it is conceivable that the minimum disk mass required to generate 
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spiral arms might exceed the maximum disk mass allowed by the rotation curve (see Quillen & 
Pickering 1997). This may provide a further test. 

It may be that spiral arms do not have an origin internal to the dynamics of the disk. 
However, it is difficult to invoke an external trigger in LSB galaxies since they are quite isolated 
systems (Bothun et al. 1993; Mo et al. 1994). If spiral structure has nothing to do with disk 
kinematics, then of course there is no test one way or the other. 

The form of the surface brightness distribution in the field (McGaugh 1996; Fig. 8 of paper I) 
may also hold a clue. Milgrom (1989) argued that the stability of disks would lead to a transition 
at the critical surface density ctq ~ a$G~ l = 88OA40pc~ 2 . Above this, the typical acceleration of 
a disk will exceed clq and one expects purely Newtonian behavior. Since bare Newtonian disks are 
subject to self destructive instabilities, they should not survive. Below the critical surface density, 
MOND lends the extra stability discussed above and a disk may assume any gq < ctq. 

There should therefore exist a sharp cut off in the surface brightness distribution at 
Eq* = ctq/Y*. Such a sharp cut off is observed, though its precise position is rather uncertain 
(McGaugh 1996). For Y* = 2, one expects /Uo* = 20.4, consistent with the bright edge of the 
surface brightness distribution determined by de Jong (1996). 

The distinction between giant ellipticals and spirals may be dictated by the critical ctq. Gas 
dissipation tends naturally to form disk systems. However, these are only stable if do < (Jq. An 
elliptical galaxy may result when, whether through initial conditions or subsequent mergers, a 
system exceeds the critical surface density. It will not be stable as a disk, but since there is no 
objection to three dimensional, pressure supported, Newtonian systems, an elliptical galaxy seems 
a natural result. The collapse of rotational support presumably has drastic consequences for the 
gaseous content of such a system, so one might expect a starburst to consume the gas roughly 
coeval with the instability event. This would leave the signature of a population dating to said 
event, so that at present elliptical galaxies would appear gas poor and old. 

This same process could also occur with conventional dynamics and dark matter halos. Should 
self-gravity become sufficiently great in high surface density disks, the halos will no longer suffice 
to stabilize them. However, the dark matter picture offers no reason why this should happen at 
the particular scale clqG~ 1 natural to MOND. 

The thickness and velocity dispersion of disks may provide a further test. In the Newtonian 
case the dark mass is arranged in a halo, so the thickness of the disk is determined solely by 
the mass of the stars and the usual conventional dynamical equations. For a vertical density 
distribution 
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(Spitzer 1942), 
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This would lead one to expect the disks of LSB galaxies either to be relatively thick or to have very 
low vertical velocity dispersions since the low surface mass density disks have little self-gravity. 
In contrast, MOND increases the binding force over the Newtonian prediction in a way which 
increases with decreasing surface density. The modified version of the expression for disk thickness 
(Milgrom 1983b) becomes 



where fi( x ) is the interpolation function of MOND. Disks in the far MOND regime should be 
thinner than the Newtonian equivalent by the factor fj,(x) « x = a/ao for a given velocity 
dispersion. For a given disk thickness, MOND disks can support a velocity dispersion a factor of 
\J fi(x) higher. 

This could provide a strong test. Low surface brightness disks appear to be quite thin 
(Dalcanton k. Schectman 1996; Kudrya et al. 1994). If they are Newtonian they must have quite 
small velocity dispersions. MOND disks of the same thickness would have distinctly higher q z . 

Moreover, there comes a surface density where Newtonian disks cease to be disks at all. To 
illustrate this, consider the disk thickness resulting from equations (27) and (28) with velocity 
dispersions plausible for the central regions of low surface brightness disks. The only difference 
between these two equations is the MOND interpolation function fi(x). For plotting convenience 
we adopt fi(x) = x/s/l + x 2 (Milgrom 1983a). The interesting effects occur in the asymptotic 
regime n(x) ~ x where the assumed form of fi(x) is irrelevant. Assuming the central regions of 
the disk can be approximated as a plane parallel slab where V 2 /R <C s%/zq so that the vertical 



restoring force dominates the acceleration allows us to approximate x as x « yao/a^. This is 
useful for illustrating the dependence of the disk thickness on the central surface brightnesses and 
velocity dispersion without stipulating a specific form of V{R). In general, the precise MOND 
prediction depends on the total acceleration, and radial, tangential, and vertical components can 
all contribute. The approximation made here using only the vertical component is adequate for 
illustrating the relevant effects near the centers of disks. The most interesting effect is the behavior 
behavior of purely Newtonian disks of very low surface mass density. 

The disk thickness is plotted as a function of disk central mass surface density in Fig. 9 for 
plausible values of the velocity dispersion and disk scale length. Note that as the central mass 
surface density decreases, there comes a point (which also depends on ? 2 ) where Newtonian disks 
rapidly become intolerably thick. Unless q z is quite low for very LSB disks, these objects should 
not be disks at all. Disks remain reasonably thin in MOND because the restoring force of the disk 
is larger. 

Consider the actual numbers for the case illustrated in Fig. 9. A normal sized (h = 3 kpc) 
disk with a central velocity dispersion q z = 20 kms -1 is respectably thin in both Newtonian and 
MOND cases for <7o > 8OA40pc~ 2 (fio ~ 23 for T* = 2). Below this surface density, the two 
predictions diverge. MOND disks remain credibly thin (zo/h < 1/4) down to do ~ 2A4qpc -2 
(Ho ~ 27). In the Newtonian case, zo/h = 1/4 occurs at ctq ~ 30A4qpc -2 (/xq ~ 24) and by 
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cro = lO.A/f0pc~ 2 (/xo ~ 25) the object ceases to be a disk at all, with zq /i. In order to have 
such low surface brightness disks (which do exist) the central velocity dispersion must be very low: 
<j 2 < 10 kms" 1 is required to postpone the Newtonian divergence to the regime /j,q > 25 not yet 
observed. Global stability is not the only problem for purely Newtonian disks: it is also difficult 
to explain the existence of cold, thin disks of low surface brightness with a purely Newtonian force 
law. 

Unfortunately, no stellar velocity dispersion data exist for LSB galaxies. This would be an 
extremely difficult observation, but would provide a powerful test. Velocity dispersion data do 
reach to large radii in some high surface brightness disks (Oiling 1996a). If interpreted in terms of 
dark matter, the MOND signature in such data will be a requirement either for a large amount 
of disk dark matter or a flattened halo. Data which reach far enough to imply a large mass 
discrepancy should make it necessary to put a lot of dark mass in a distribution close to that of 
the disk (e.g., Oiling 1996b). In very LSB disks, there may come a point where the dark matter 
required to bring this about would exceed that allowed by the maximum disk solution. 

A related observation is the Oort discrepancy in the Milky Way. Kuijken & Gilmore (1989) 
test MOND in this context, and find that it tends to over-correct somewhat. However, they used 
a value of ao nearly four times larger than currently measured. This will cause a mass discrepancy 
to be implied by the MOND equations before it actually should be, leading to an apparent 
over-correction. 

3.4. Clusters of Galaxies 

Perhaps the strongest observational argument against MOND at the time of its introduction 
was that while it worked in galaxies, it failed in clusters of galaxies. The apparent mass 
discrepancies Mn/M* ^ 100 were too large to be explained by the typical accelerations, 
<W( a ) ~ 10- Even with MOND, substantial amounts of dark matter seemed to be required, an 
unacceptable situation. 

The basic picture which held for many years is that the amount of dark matter increased 
with increasing scale. Mass discrepancies in galaxies were a factor of ten, those in clusters factors 
of hundreds, and even more was required to close the universe. This picture of ever more dark 
matter on ever larger scales has changed dramatically in recent years. X-ray observations of rich 
clusters of galaxies have shown that much or even most of the baryonic matter is in the form of 
hot gas (e.g., David et al. 1990; White & Fabian 1995). Rather than having more dark matter 
than individual galaxies, clusters are in fact more baryon dominated (paper I). There were in 
effect two missing mass problems in clusters: the usual dynamical mass discrepancy and the fact 
that many of the baryons are in the form of hot gas rather than stars. 

The need for an additional mass component which seemed so unacceptable for MOND would 
now appear to be a successful prediction. It is thus important to reanalyze MOND in light of 
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this new knowledge. The & White (1988) did this for the case of the Coma cluster, and found 
that MOND still seemed to be off by a factor of ~ 2. Recall, however, the many uncertainties 
associated with tests in quasi-spherical systems discussed in the case of dwarf Spheroidals. Given 
the uncertainty in the underlying assumptions of sphericity and virial equilibrium upon which the 
analysis is based, is a factor of two a problem or a success? 

Andernach & Tago (1997) have compiled a good deal of cluster redshift data, and quote a 
median cluster velocity dispersion of 695 kms -1 . For a spherical virialized system, this corresponds 
to a MOND mass of 3 x 10 14 A4q (equation 22). The typical hot gas mass of clusters in the 
compilation of White & Fabian (1995) is ~ 10 14 A4q. So even a crude calculation comes pretty 
close, especially if there is a comparable amount of mass in stars. 

Sanders (1994) re-addressed the problem in greater detail. For 16 clusters, he shows that the 
MOND mass is strongly correlated with the X-ray gas mass. The MOND masses are a bit greater 
than the gas mass (his Fig. 3) indicating some additional mass in stars. 

A more thorough investigation of this particular test would require a great deal more data 
than are available to us: combined X-ray and optical observations giving good estimates of the 
mass in each component and accurate velocity dispersions. Obtaining the full optical luminosity 
of a cluster is far from trivial (Impey et al. 1988), and X-ray observations still lack the spatial 
resolution to guarantee the validity of the assumption of isothermality. Obtaining a reliable 
velocity dispersion from a sample of galaxies that are truly virialized is also challenging. Interlopers 
and substructure could play a strong role is distorting the mass-indicative velocity dispersion. A 
40% overestimate of q would imply a MOND mass a factor of 4 too great. This would appear 
disastrous if this uncertainty is not considered and properly propagated. It is important to test 
MOND in clusters as adequate data become available, keeping in mind the many uncertainties. 

3.5. Gravitational Lensing 

Another indication of a mass discrepancy independent of kinematic data is gravitational 
lensing, both in individual galaxies (e.g., Kochanek 1995) and clusters (e.g., Tyson et al. 1990). 
Since there is not yet a relativistic extension of MOND, there is no clear prediction for gravitational 
lensing. That lensing occurs must ultimately be explained, but at present it provides no test of 
the validity of the MOND force law. 

Some progress has nevertheless been made. If MOND is interpreted as an alteration of the 
law of inertia, it is fairly successful at explaining weak lensing in clusters. However, it does then 
predict that additional, as yet undetected baryonic mass resides in the cores of a few of the 
most X-ray rich clusters (Milgrom 1996). This is apparently the only place where MOND does 
not remedy the mass discrepancy problem in the sense that significant additional mass remains 
hidden. This is a very limited missing mass problem, restricted to the central Newtonian cores of 
some rich clusters. These are generally cooling flow clusters, so at least some additional baryonic 
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mass is already inferred to reside there. 

Gravitational lensing observations have allowed entire classes of theories to be ruled out 
as the possible basis of combining relativity and MOND (Bekenstein & Sanders 1994). This 
does not mean none can exist. One attempt at such a theory (Sanders 1997) has an interesting 
consequence. In this case, the mass distribution that would be inferred when interpreted in terms 
of dark matter is the same for both lensing and kinematical observations. Hence, observational 
agreement between such observations does not uniquely require a dark matter interpretation. 

3.6. Large Scale Structure 

Observations of the motions of galaxies now extend over very large scales. These do not 
provide strong tests of MOND, since we must make crude assumptions about the mass distribution. 
We can, however, address the consistency of MOND with current observations on a qualitative 
level. 

That there are large (;> 30 Mpc) voids and filaments in the large scale galaxy distribution is 
now regularly reproduced by simulations of structure formation. It is worth recalling that these 
features initially came as a great surprise (de Lapparent et al. 1986): the universe was suppose to 
be homogeneous on these large scales. That these large, sharp features exist makes a reasonable 
amount of sense in the presence of a force law which is effectively enhanced on large scales, as does 
the occasional inference of an excessive mass concentration like the Great Attractor (Lynden-Bell 
et al. 1988). 

Where dynamical measurements exist, these large scale structures can be used as crude tests 
of MOND. Milgrom (1997) presents an analysis of the Perseus-Pisces filament, and derives T ~ 10. 
This is fairly reasonable considering the crude assumptions that must be made (i.e., that the 
structure is a virialized linear feature) and that it may well contain substantial amounts of gas. 

The expansion of voids could provide another test. As underdense regions, voids must expand 
in the conventional picture. This need not be the case in MOND, for which the effective potential 
is logarithmic on large scales. A spherical shell mass distribution will eventually turn around 
and collapse regardless of its interior density (Felten 1984), so one might expect to see some 
voids expanding and others contracting. In the one case where careful distance as well as redshift 
measurements have been made (Bothun et al. 1992), no expansion is detected (^ 5% of the void 
diameter in kms -1 ). 

On the largest scales, MOND does require a universe composed entirely of baryons. It 
may therefore seem troubling that dynamical estimates of the mass density are persistently 
around 0, ~ 0.3, much higher than the baryon density allowed by primordial nucleosynthesis, 
0.01 < < 0.03. These estimates of Q are of course based on the usual Newtonian equations, 
and MOND will require much less mass. As with all mass discrepancies, the amount by which Q 
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is overestimated depends on the typical acceleration scale probed: 

^MOND = — fi- (29) 

a 

It is difficult to estimate (a) here, but as an example we make the usual assumption of 
homogeneity so that the density field may be approximated as a constant. The Poisson equation 
with the usual Q = 8ttGp/3Hq then gives 

(a) « \V<p\ = ^HfaR, (30) 

where R is the scale which the observations cover. Davis et al. (1996) give a value ri 6 6 = 0.6 
from data which extend reliably out to R « 5000 kms -1 . There should be no mass biasing in 
MOND (though different populations of galaxies may be biased relative to one another) so we 
assume 6=1. This gives a conventional = 0.43. Taking these numbers at face value leads to 
(a) « 0.026As~ 2 and 

^MOND ~ 0.01. (31) 

Given the nature of the data and the necessary assumptions, this is probably uncertain by at 
least a factor of a few. Nevertheless, it is striking that 

^mond ~ £l b . (32) 

MOND appears to adequately address the dynamical mass discrepancy problem on even the 
largest scale. 



4. Discussion 

We have taken care to review previous analyses which have found fault with MOND. There 
is no evidence on any scale which clearly contradicts MOND. Some data which are cited as 
contradicting MOND actually appear to support it (e.g., that for dwarf Spheroidals). It should be 
noted that not all previous independent analyses of MOND have been negative. Begeman et al. 
(1991), Morishima & Saio (1995), and Sanders (1996) all report positive tests of the theory. The 
new data for LSB galaxies we have collected obey all the predictions made by Milgrom (1983b). 

The observational tests of MOND we have discussed are summarized in Table 3. That 
LSB galaxies fall on the Tully-Fisher relation is a strong prediction of MOND (§2.3). The T-S 
conspiracy which occurs when interpreting the LSB galaxy Tully-Fisher relation with dark matter 
(Zwaan et al. 1995; paper I) can be derived from MOND (§2.4). Stellar mass-to-light ratios which 
agree well with what is expected for stars can be computed directly from the observations with 
MOND (§2.5). The mass surface density implied by the kinematic measure V^/h is strongly 
correlated with the disk surface brightness with the slope expected in MOND (§2.6). The 
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radial variation of the dynamical mass-to-light ratio computed conventionally behaves in a 
manner consistent with the predictions of Milgrom (1983b). Similarly, the radii of transition to 
apparent dark matter domination and the typical accelerations observed in disks vary with surface 
brightness as expected in MOND (§2.7). The systematic dependence of the shape of rotation 
curves on surface brightness is also predicted by MOND (§2.8). Taken in sum, the data are well fit 
by MOND (§2.9; paper III). Indeed, all of the empirical facts we identified in paper I describing 
the systematic properties of the rotation curves of disks as a function of surface brightness were 
anticipated by Milgrom (1983b). 

MOND also survives tests in systems other than disks. Dwarf Spheroidal galaxies are an 
important example. Once thought to fail there (Gerhard & Spergel 1992), MOND now appears to 
do well with improved data (§3.1). On the other hand, giant elliptical galaxies provide no useful 
test of MOND since the accelerations they experience do not probe the MOND regime, at least 
not with current observations (§3.2). The stability properties of disks appear consistent with 
MOND, and the velocity dispersions of thin LSB disks could provide a very strong test (§3.3). 
Galaxy clusters provide another important test. Originally seeming to require additional dark 
matter in clusters (Milgrom 1983c), the detection of large amounts of hot X-ray emitting gas in 
clusters generally brings these into consistency with MOND (§3.4; Sanders 1994). A definitive 
prediction for gravitational lensing requires a relativistic generalization of MOND which does not 
yet exist (§3.5). Some types of theories can be ruled out on this basis (Bekenstein & Sanders 1994) 
while others remain possible (Sanders 1997). On the largest scales, MOND does as well as can be 
expected given the applicability of the available data (§3.6). A very low density (f2 ~ 0.01), purely 
baryonic universe is roughly consistent with the dynamical data which constrain £1. 

Empirically, MOND is the effective force law in disk galaxies. It appears that this may also 
be the case in other systems. The reason for this phenomenology needs to be understood. 

There data allow two possible interpretations. Either 

1. MOND is correct, or 

2. Dark matter mimics the behavior of MOND, at least in disks. 

The second possibility implies a unique and powerful coupling between dark and luminous matter. 
It is possible to write down an equation which directly links the dark matter dominated dynamics 
to the detailed distribution of the luminous matter. This provides a new observational test of 
theories of disk galaxy formation within the standard dark matter paradigm. Since MOND 
always fits disk galaxy rotation curves, it must be possible to take the luminous mass distribution 
predicted for any given disk by a dark matter galaxy formation theory, apply the MOND procedure 
to the luminous mass only, and thereby obtain the correct rotation curve. If this can not be done, 
the theory has failed to produce a realistic disk. 

There are no clear empirical objections to the first possibility. Milgrom (1983b) did accurately 
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predict numerous aspects of the kinematical properties of LSB galaxies. This seems unlikely to 
have occurred by accident, so the possibility that MOND is correct should be considered seriously 

We are grateful to Moti Milgrom, Chris Mihos, Vera Rubin, Bob Sanders, and the referee for 
close reading of this manuscript and numerous helpful conversations. 
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Table 1. MOND T* 



Galaxy 


I* 




F563- 


-1 


2.6 


0.8 


F563-V2 


1.4 


1.6 


F568- 


-1 


1.2 


0.7 


F568- 


-3 


1.8 


1.0 


F568- 


-VI 


3.7 


1.3 


F571- 


-VI 


-0.2 a 


0.4 


F574- 


-1 


0.01 b 


0.6 


F583- 


-1 


0.9 


1.3 


F583- 


-4 


0.2 


0.5 


UGC 


128 


1.0 


0.5 


UGC 


6614 


2.7 


0.7 


DDO 


154 


~ C 


1.2 


DDO 


168 


0.8 


0.9 


NGC 


55 


0.3 


1.3 


NGC 


247 


2.2 


1.2 


NGC 


300 


1.6 


1.2 


NGC 


801 


1.2 


1.1 


NGC 


1560 


2.0 


1.3 


NGC 


2403 


1.5 


2.4 


NGC 


2841 


~ 7 C 


6.1 


NGC 


2903 


4.8 


5.5 


NGC 


2998 


0.8 


2.3 


NGC 


3109 


0.4 


0.8 


NGC 


3198 


2.5 


2.6 


NGC 


5033 


5.4 


2.3 


NGC 


5533 


4.1 


1.8 


NGC 


5585 


0.9 


1.6 


NGC 


6503 


1.7 


2.3 


NGC 


6674 


2.9 


2.3 


NGC 


7331 


2.6 


3.5 


UGC 


2259 


2.7 


1.7 


UGC 


2885 


1.6 


1.8 



a ±2.5 
b ±1.4 

c distance sensitive 
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Table 2. Dwarf Spheroidals 



Galaxy 


r* 


V 


Carina 


5.4 


0.5 


Draco 


10.8 


0.9 


Fornax 


0.4 


0.6 


Leo I 


2.1 


3.8 


Leo II 


2.5 


1.4 


Sculptor 


1.2 


0.6 


Sextans 


2.9 


0.3 


Ursa Minor 


16.9 a 


0.5 



a T* = 1012 for an 
alternative measurement 
of?. 
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Table 3. Tests of Predictions 



Observational Test 


MOND 


LSBG Tully-Fisher Relation 


V 


T-Y, Relation 


V 


Stellar Mass-to-Light Ratios 


V 


Mass Surface Densities 


V 


Conventional T a (R) 


V 


Transition Radii 


V 


Characteristic Accelerations 


V 


Rotation Curve Shapes 


V 


Rotation Curve Rate of Rise 


V 


Rotation Curve Fits 


V 


Disk Stability 


? 


Dwarf Spheroidal Galaxies 


V 


Giant Elliptical Galaxies 


NT 


Galaxy Clusters 


? 


Gravitational Lensing 


NP 


Large Scale Structure 


? 


q = n b ? 


? 



Note. - 

y/ = prediction confirmed 
X = prediction falsified 
? = remains uncertain 
NP = no prediction 
NT = no test 
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Fig. 1. — The stellar mass-to-light ratios of spiral galaxies determined from MOND, plotted against 
(a) absolute magnitude, (b) central surface brightness, and (c) scale length. The median value is 
Y* = 1.6. This Figure is analogous to Fig. 4 of paper I. As was done there, error bars in Y* 
have been computed assuming a nominal inclination uncertainty of 3 degrees. The errors are larger 
than those in Fig. 4 of paper I because MOND masses are proportional to the fourth power of 
the velocity, and so have a factor sin _4 (i) contributing to the error instead of just sin~ 2 (i) in the 
conventional case. 

Fig. 2. — The mass surface density £ indicated by MOND plotted against surface brightness. The 
two are strongly correlated, as expected in the absence of dark matter. The line is not a fit to the 
data. It illustrates the slope expected from the form of the MOND force law. The normalization 
is also predicted modulo the baryonic mass-to-light ratio Y^ (equation 20). The line is drawn for 
Y 5 = 4. 

Fig. 3. — The accumulated conventional dynamical mass-to-light ratio (Y D oc V 2 R/L) as a function 
of radius. The plot is divided into bins of different central surface brightness as indicated in each 
panel. All available data (Tables 1 and 2 of paper I) have been used. In MOND, one expects the 
mass discrepancy indicated by Y Q to be larger and to set in at smaller radii in galaxies of lower 
surface brightness. 

Fig. 4. — The rotation curves plotted in terms of the requisite centripetal acceleration a = V 2 (R)/R. 
The plot is divided into different bins in central surface brightness, and the value of ao is marked 
by the dashed line. In MOND, one expects the acceleration to decline with surface brightness. 

Fig. 5. — The rotation curves plotted logarithmically to illustrate their shapes. This test is not 
inclination dependent. In MOND, one expects high surface brightness galaxies to have rapidly 
rising rotation curves which fall gradually to the asymptotically flat value. Low surface brightness 
galaxies should have rotation curves which rise slowly to the asymptotically flat value. 

Fig. 6. — An example of the MOND fits to the rotation curves of low surface brightness galaxies 
from paper III. The data points are for the galaxy F563-1. The dashed line shows the Newtonian 
rotation curve of the stellar disk and the dotted line that of the gas. The solid line is the resulting 
MOND fit to the entire rotation curve. This follows directly when the MOND force law is applied 
to the observed luminous mass distribution. 

Fig. 7. — The residuals (in percent) of all the MOND fits to the rotation curves of LSB galaxies. The 
difference between the observed velocity and that predicted by MOND from the observed luminous 
mass distribution (AV = V t, s — Vmond) is shown as a function of the critical parameter x = a/ao 
where in practice a = V 2 /R. Each point represents one measured point from the rotation curve 
fits described in paper III. The residuals for all points from all LSB galaxies are shown together. 
This represents 15 new LSB galaxy fits with nearly 100 total independent measured points. 

Fig. 8. — The MOND mass-to-light ratios inferred for Local Group dwarf Spheroidal galaxies. 



-35- 



Triangles are the determinations of Gerhard (1994). Solid triangles are the values he finds for the 
isolated case appropriate if 77 > 1 (equation 24) , and open triangles are his quasi-Newtonian values 
which are more appropriate when the external field is dominant (77 < 1). Lines capped by crosses 
illustrate the range of allowable values determined by Milgrom (1995). The results of Gerhard and 
of Milgrom have been scaled to the value of ao adopted here which increases their Y* by a factor of 
1.6. Solid circles are our own determinations based on the data compiled by Mateo (1996). Dashed 
lines delimit the most plausible range, 1 < Y* < 6. Most galaxies fall in this range according to all 
three independent determinations. In only one case (Ursa Minor) is there a marginally significant 
deviation from the most plausible range. 

Fig. 9. — The thickness zo/h expected for disks of various central surface densities do- Shown along 
the top axis is the equivalent .B-band central surface brightness /Jo for Y* = 2. Parameters chosen 
for illustration are noted in the Figure (a typical scale length h and two choices of central vertical 
velocity dispersion q z ). Other plausible values give similar results (equations 27 and 28). The 
solid lines are the Newtonian expectation and the dashed lines that of MOND. The only dif ference 
between the two cases is the MOND interpolation function //(x) = x/\/l + x 2 with x = ^(Tq/oq. 
The Newtonian and MOND cases are similar at high surface densities but differ enormously at 
low surface densities. Newtonian disks become unacceptably thick unless LSB disks are very cold 
{q z < 10 kms" 1 ). In contrast, MOND disks remain reasonably thin to quite low surface density. 
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